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Comment on "Unconventional Magnetism in a 
Nitrogen-Containing Analog of Cupric Oxide" 

In a recent Letter. Zorko et al. [l[ report on an "un- 
expected inhomogeneous magnetism" related to "a pe- 
culiar fragility" of the resonating-valence-bond (RVB) 
electronic state in the spin-i quantum magnet CuNCN. 
Here, we show that: i) the spin model of the anisotropic 
triangular lattice (ATL), as proposed in Ref. does 
not apply to CuNCN; ii) this model does not support the 
proclaimed RVB ground state in the relevant parameter 
range; iii) the magnetic state of CuNCN is not necessarily 
inhomogeneous. 

In Ref. the ATL geometry is derived from the two 
shortest Cu-Cu distances in the crystal structure (Fig. [TJ 
left panel) and does not account for the 7r-conjugatcd 
NCN groups that mediate strong long-range superex- 
change interactions. The sizable dispersion of Cu 3d 
bands along the c direction (Fig. 3 in Ref. Q) indeed im- 
plies a large electron hopping and, consequently, a strong 
antiferromagnetic (AFM) exchange J ~ 2500 K medi- 
ated by the NCN groups (Fig. [TJ right panel). The sim- 
ilar NCN-driven long-range coupling has been reported 
in closely related transition-metal carbodiimides Af NCN 
with M = Fe, Co, and Ni (3, Q. The strong coupling 
J is by far the leading energy scale of the system that 
does not simply modify the ATL physics, but essen- 
tially changes the microscopic scenario. The quasi- two- 
dimensional J\ — Jb problem in the ab plane is replaced 
by the quasi-one-dimensional (ID) lattice of spin chains 
running along the c direction (see Fig.[TJ. 

Even if the leading coupling J along the c direction 
is ignored, the physics in the ab plane does not support 
the proclaimed RVB scenario. According to the precise 
numerical study based on the coupled-cluster method 
(CCM) m, the phase space of the spin-i Heisenberg 
model on the ATL lattice features long-range-ordered 
(LRO) phases only. Other theoretical studies consis- 
tently restrict the possible parameter range of the pro- 
posed gapped RVB-type phase to 1.2 < J\/ Jb < 1.6 0, 
while Zorko et al. estimate J 1 = 2300 K, J b = 540 K, 
and Ji/Jb — 4.3, which is well inside the collinear stripe- 
ordered phase of the ATL model 0, Moreover, the 
leading coupling Ji of 2300 K is found for the Cu-N-Cu 
angle of 96.75°, where a ferromagnetic or weakly AFM in- 
teraction is expected from Goodenough-Kanamori rules. 
An RVB scenario also contradicts the computational es- 
timates of Ref. ( Ji ~ -500 K, J b ~ 10 K). 

The RVB ground state lacks static magnetic fields and 
shows short-range spin correlations only. However, both 
muon spin rotation (/iSR) and nuclear magnetic res- 
onance (NMR) experiments evidence the formation of 
static fields below 70 — 80 K. To circumvent this dis- 
crepancy, the authors of Ref. [l[ suggest the "peculiar 
fragility" of the RVB state, which is allegedly destroyed 
by the muons and by the magnetic field of 9.4 T in the 
/iSR and NMR experiments, respectively. In this situ- 



ation, only the electron spin resonance (ESR) measure- 
ment performed at a lower field of 0.1 T could access the 
true ground state, but the ESR response is at best am- 
biguous and does not show any specific footprints of the 
RVB scenario. The decrease in the ESR intensity below 
70 K can be evenly understood as the onset of the conven- 
tional long-range magnetic order that is indeed expected 
from microscopic considerations [H. 

Starting from this microscopic viewpoint, we interpret 
the static magnetic fields observed in the ^SR and NMR 
experiments as further signatures of the LRO state below 
70 K. In /iSR, the broad distribution of internal fields can 
be ascribed to several incquivalent muon stop sites, while 
the low value of the stretch exponent (inset to the right 
panel of Fig. 3 in Ref. [H) is well understood in terms 
of two distinct relaxation times Therefore, neither 
/iSR nor NMR show unambiguous signatures of an inho- 
mogeneous magnetic state. The experimental data can 
be thus reconciled in terms of the conventional LRO sce- 
nario and do not require an unnecessary assumption of 
an unconventional fragile electronic state. 

In conclusion, the RVB ground state of CuNCN is not 
supported by direct experimental evidence, and contra- 
dicts both: i) well-established numerical results for the 
ATL model [5J-17J; ii) microscopic description of the com- 
pound in terms of its electronic structure 0. Further 
work on CuNCN is presently underway to better under- 
stand experimental observations. 

Alexander A. Tsirlin and Helge Rosner 
MPI-CPfS, Dresden, Germany 
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FIG. 1. Microscopic magnetic model of CuNCN: anisotropic 
triangular lattice (ATL) proposed in Ref. [J] (left panel), and 
uniform spin chains running along the c direction, according 
to band structure calculations @] (right panel). The leading 
exchange coupling is J ~ 2500 K compared to Ji ~ —500 K 
and J b < 2 K. 
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Addendum to the Comment 

The magnetic state of CuNCN below 70 K has been 
a matter of debate in the recent literature [IHll- Here, 
we provide extra comments on publications preceding the 
Letter by Zorko et al. [H , and discuss their results in the 
context of the microscopic magnetic model from Ref. [l| . 
We focus on two main problems: i) the microscopic mag- 
netic model of CuNCN; ii) experimental observation of 
the potential long-range-ordered (LRO) magnetic state 
in this compound. The following comments extend our 
published work [l[ , and clarify its relation to experimen- 
tal results presented in Refs. 0, Q. 

Reliable information on the electronic structure is 
a necessary prerequisite of the microscopic magnetic 
model. The electronic structure of CuNCN has been 
first considered by Liu et al. 0, who performed density 
functional theory (DFT)+E/ calculations and arrived at 
the anisotropic triangular lattice (ATL) spin model (see 
Fig. 1 of the Comment). Their results apparently con- 
tradict band dispersions shown in Fig. [21 because the 
strong antiferromagnetic (AFM) coupling Jf, is inconsis- 
tent with the nearly flat Cu bands along the b* direction 
of the reciprocal space. By contrast, the large band dis- 
persion along c should lead to a strong AFM coupling, 
which is, however, not taken into account in the ATL 
model of Ref. 0. The ATL model was indeed disproved, 
but - surprisingly - for a totally different reason, when 
Xiang et al. [3| claimed a non-magnetic state of Cu +2 
based on new band structure calculations. 

The rather unconventional conclusion on the non- 
magnetic Cu +2 in an insulating compound contradicts 
the ESR response observed by Zorko et al. W. This Let- 
ter further refers to the ATL model of Ref. |2} , although 
its results were disproved in Ref. [3 by the same group. 
To the best of our knowledge, the relationship between 
the mutually contradicting results of Refs. 0, Q and the 
model assumptions of Ref. [I| has not been clarified. 

Once all relevant interactions are taken into account, 
DFT-based calculations provide a consistent microscopic 
picture that is robust with respect to subtle technical 
details, such as the treatment of strong electronic corre- 
lations on the Cu +2 site. The conclusion on the strong 
AFM coupling J mediated by the NCN groups stems 
from the large band dispersion along the Y — Z direction 
of the reciprocal space (see Fig. [2]). By contrast, the dis- 
persion along r — Y is vanishingly small, thus suggesting 
the negligible role of J& that, according to Refs. |1, E[, is 
an integral part of the anisotropic triangular spin lattice. 
These results are robust with respect to a specific choice 
of the exchange-correlation potential, and have been fur- 
ther confirmed by DFT+C/ calculations The micro- 
scopic analysis also puts forward the close similarity be- 
tween CuNCN and Cu +2 oxides, because the relevant lig- 
and states of the NCN unit are formed by 2p orbitals of 
nitrogen, similar to oxygen 2p orbitals in cuprates. Con- 
sidering the long-standing experience of successful ap- 
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FIG. 2. Left panel: band structure of CuNCN calculated 
in local density approximation (LDA) and the fit with the 
tight-binding (TB) model shown in thin orange and thick 
green lines, respectively. Circles denote the contribution of 
the half-filled Cu d x 2_ y 2 orbital. The large dispersion along 
Y — Z identifies the leading AFM coupling J along c, while flat 
bands along Y — Y clearly show the very weak AFM coupling 
Jt- Finally, the complex dispersion along Z — T suggest the 
competition of the nearest-neighbor (Ji) and next-nearest- 
neighbor (J2) couplings along the structural CuN2 chains (a 
direction). Note that in this analysis we start from the un- 
corrected (metallic) LDA band structure and obtain the in- 
sulating solution after mapping the LDA band structure onto 
an effective Hubbard model. 



plications of DFT and DFT+C/ to Cu+ 2 oxides @, the 
microscopic scenario for CuNCN is reliable. In a car- 
bodiimide system, an unexpected failure of DFT-based 
methods that have proved their efficiency for oxides can 
be ruled out. Therefore, the leading AFM coupling J has 
to be taken into consideration, and the ATL spin model 
of CuNCN should be replaced by the microscopic-based 
one (Fig. 1 of the Comment). 

Concerning experimental results, the onset of static 
magnetic fields at 70 K is neither accompanied by an 
anomaly in the specific heat, nor observed by neutron 
scattering. Zorko et al. 0| consider these experimental 
facts as additional arguments favoring the "unexpected 
inhomogeneous magnetic state" rather than the conven- 
tional LRO in CuNCN. The onset of the magnetic or- 
der without the corresponding specific-heat anomaly and 
neutron signal is surprising indeed, although not uncom- 
mon in low-dimensional spin-i quantum magnets with 
low ordered magnetic moments, caused by strong quan- 
tum fluctuations, and diminutively small magnetic en- 
tropy at the Neel temperature Tjv- 

To better understand the experimental results, we per- 
formed further calculations for the microscopic magnetic 
model of CuNCN established in Ref. [f[. Our quan- 
tum Monte-Carlo simulations @ reveal T)v/J — 0.12 and 
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FIG. 3. Spin lattice of CuNCN according to the compu- 
tational results of Ref. J ~ 2500 K, Ji ~ -500 K, 
J 2 ~ 100 K, Jac ~ 70 K, and J b < 2 K. Open and filled 
circles denote up and down spins in the ground state spin 
configuration derived from CCM calculations. The ordered 
magnetic moment is about 0.4 \ib- 
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FIG. 4. Magnetic specific heat of CuNCN calculated with 
quantum Monte-Carlo simulations, and the respective mag- 
netic entropy shown as an inset. Note the absence of any visi- 
ble anomaly at the magnetic ordering transition Tjv/J ~ 0.12 
due to the very low entropy available at Tjv. 



confirm the absence of any visible anomaly in the mag- 
netic specific heat (Fig. @|. Experimentally, the mag- 
netic contribution to the specific heat is superimposed 
onto the lattice contribution (about 20 J mol -1 K _1 at 
Tn = 70 K) that further impedes the experimental obser- 
vation of the magnetic transition. The small transition 
anomaly originates from the low magnetic entropy (be- 
low 3 % of the full magnetic entropy R In 2) available at 
Ttv. Similar results for the missing specific-heat anomaly 
at Tjv have been recently obtained for a family of quasi- 
two-dimensional spin-i quantum magnets 0. 

The absence of the magnetic neutron scattering can be 
ascribed to a spin-liquid ground state, which is devoid of 
the long-range magnetic order, or to the LRO state with 
the small magnetic moment that lies below the sensitivity 
threshold of the neutron-scattering experiment. Apply- 
ing the coupled-cluster method to our DFT-based micro- 
scopic model, we establish the stripe AFM ordering pat- 
tern (Fig. [3]) and estimate the ordered magnetic moment 
of about 0.4 hb in CuNCN. An experimental observa- 
tion of such a low magnetic moment remains challenging 



and requires instruments with excellent sensitivity and, 
preferably, single crystals, while in Refs. Q and [3[ none 
of these conditions were fulfilled. We believe that the 
interpretation of the neutron data is premature, because 
neither Liu et al. Q nor Xiang et al. Q compared the 
sensitivity of their neutron measurements to the magni- 
tude of the anticipated magnetic scattering in CuNCN. 
The systems showing the LRO magnetic state and lack- 
ing any magnetic scattering observable in a powder neu- 
tron experiment include, for instance, Sr2Cu03 Q and 
PbV03 0. In quantum magnets, this situation is not 
uncommon, and thus the reported neutron experiments 
alone can not be a decisive test for the formation of the 
LRO magnetic state. 

We would like to acknowledge Ronald Zinke and Jo- 
hannes Richtcr for providing the CCM results. 

Alexander A. Tsirlin and Hclge Rosncr 
MPI-CPfS, Dresden, Germany 
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